AD-A170  488 


>■  i  ^  ■  U  ■  1 1  V  9  \,  ■  IJ  ^  I  ■!<  V  n  <.^1^5  tf  ■ 'J 


Department  of  Mechanical 
and  Industrial  Engineering 

University  of  Illinois  at  Urbana-Champaign 
Urbana,  IL  61801 


EFFECT  OF  CONVECTION  ON  WELD  POOL  SHAPE  AND  MICROSTRUCTURE 


J.  Mazumder,  M.  M.  Chen,  C.  L.  Chan  and 
R.  Zehr 


EFFECT  OF  CONVECTION  ON  WELD  POOL  SHAPE  AND  MICROSTRUCTURE 


Report  to  U.S.  Navy 
(Contract  No.  N00014-84-K-0315) 


by 


J,  Mazumder,  M.  M.  Chen,  C.  L.  Chan,  and  R.  Zehr 
Department  of  Mechanical  and  Industrial  Engineering 
University  of  Illinois  at  Urbana-Champaign 
1206  West  Green  Street 
Urbana,  IL  61801 


i  ;  . ’ .  -.1 

i  . 

July  1986 


This  d<- 


)s»«Q  l3i/pCOT«i|  I 

c  ;  3.  IM  I 


I 


n 


w. 


V 


y  ABSTRACT 

The  role  of  convection  In  weld  pools  is  of  great  interest  due  to  its 
effect  on  weld  pool  shapes  and  weld  defects.  A  quantitative  understanding  of 
convection  in  the  weld  pool  is  being  developed,  at  the  University  of  Illinois 
at  Urbana-Champaign  under  a  program  supported  by  the  Office  of  Naval  Research 
(Contract  No.  N00014-84-K-0315) .  This  report  describes  the  progress  for  the 
past  two  years  effort  during  the  period  of  April  15,  1984  to  April  14,  1986. 

'>In  this  part  of  the  effort  major  emphasis  has  been  to  model  surface 
tension  driven  flow.  The  Lorentz  force  effect  is  incorporated  towards  the  end 
of  the  two  year  effort.  The  significant  achievements  of  the  program  are: 

a)  Stagnation  flow  analysis  for  center  of  the  weldpool  which  provided 
algebraic  expressions  for  temperature  and  velocity.  This  solution 
clarifies  the  scaling  problems. 

b)  Three-dimensional  axisymmetric  model  for  surface  tension  driven  flow 
with  flat  surface  assumption  which  predicted  the  pool  shape  for 
stationary  heat  source  such  as  in  the  case  of  spot  welds. 

c)  Three-dimensional  perturbation  model  for  surface  tension  driven  flow 
with  flat  surface  which  predicted  velocity  field,  temperature  field, 
effect  of  trace  elements  and  self-consistent  prediction  of  the  pool 
shape  and  cooling  rate. 

d)  Three-dimensional  model  with  flat  and  free  surface  which  predicted 
same  as  above  with  more  realistic  boundary  condition.  This  model  can 
also  predict  effect  of  Lorentz  Force. 

A  series  of  controlled  experiments  has  been  carried  out  to  verify  pool 
shape  and  liquid  velocity.  This  verified  the  trend  predicted  by  the  model. 
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1.  INTRODUCTION 

When  a  concentrated  heat  source,  such  as  laser,  arc,  or  electron  beam,  is 
directed  at  the  surface  of  a  moving  substrate,  several  phenomena  occur 
simultaneously.  Heat  is  being  reflected,  absorbed,  and  re-radiated  on  the 
surface.  The  absorbed  heat  develops  a  molten  region.  Inevitably,  there 
exists  a  fluid  flow  within  this  molten  region.  In  addition  to  the  convection 
due  to  the  motion  of  the  substrate,  there  is  also  the  convection  due  to  the 
fluid  flow  within  the  molten  pool.  Depending  on  the  type  of  heat  source,  the 
fluid  flow  may  be  caused  by  the  surface  tension  gradient  or  the  Lorentz  force 
or  both.  Buoyancy  force  can  also  play  a  role  in  establishing  the  fluid  flow 
in  the  case  of  casting.  When  the  dimension  of  the  molten  region  is  small  as 
in  the  case  of  welding  and  laser  materials  processing,  buoyancy  force  becomes 
unimportant.  Beyond  the  molten  region,  the  absorbed  heat  is  conducted  into 
the  substrate  and  transported  away  by  the  motion  of  the  substrate. 

A  quantitative  understanding  of  convection  in  the  weld  pool  is  being 
developed  for  the  last  two  years  at  UIUC  under  a  research  program  sponsored  by 
ONR.  In  the  last  two  years,  the  significant  achievements  of  the  program  are: 

a.  Stagnation  flow  analysis  for  center  of  the  weld  pool  which  provided 
algebraic  expressions  for  temperature  and  velocity  (47,48).  This 
solution  clarifies  all  the  scalings  of  the  problem,  i.e.  the  effects 
of  process  parameters  on  the  magnitude  of  velocity  fields,  temper¬ 
ature,  and  the  extent  of  the  viscous  effect  can  be  better  understood. 

b.  Three-dimensional  axi symmetric  model  for  surface  tension  driven  flow 
with  flat  surface  assumption  which  provided  self-consistent  pre¬ 
diction  of  the  pool  shape  for  stationary  heat  source;  velocity  and 
temperature  fields  are  also  obtained  [49). 


c.  Three-dimensional  perturbation  model  for  surface  tension  driven  flow 

with  flat  surface  [50,51]: 

•  Velocity  field  is  obtained:  Qualitative  understanding  of  the 
mixing  mechanism  has  been  achieved  by  particle  tracing 

•  Temperature  field  is  obtained  which  provided  self-consistent 
prediction  of  the  pool  shape  and  prediction  of  the  cooling  rate 

•  Effect  of  trace  element  (positive  surface  tension  coefficient) 
on  velocity  field  and  pool  shape  is  also  examined. 

NOTE:  Due  to  the  perturbation  technique  used,  this  model  is  valid 

only  for  moderate  and  low  scanning  speed. 

d.  Full  three-dimensional  model  for  surface  tension  driven  flow  with 

flat  surface  and  free  surface  [51,46]: 

Numerical  Solution 

•  Results  are  comparable  to  three-dimensional  perturbation  model; 
however,  this  model  is  capable  of  handling  high  scanning  speed 
as  well  as  low  scanning  speed 

•  Predictions  of  pool  shape,  cooling  rate  and  free  surface  shape 
These  analyses  improved  the  physical  understanding  of  the  process  and 

certain  general  features  of  the  flow  (see  Fig.  1.1).  Since  these  flow 
processes  appear  to  have  singularities  (Region  3,  Fig.  1.1)  and  certain 
limiting  behaviors,  much  work  must  be  done  to  fully  understand  the  convection 
in  the  weld.  Free  surface  and  effect  of  Lorentz  forces  should  also  be 
considered.  This  report  discusses  the  effect  of  surface  tension  driven  flow 
in  detail  and  provides  initial  data  for  free  surface  and  effect  of  Lorentz 
force  acting  simultaneously  with  surface  tension  driven  flow. 


3.  Cooled  Corner  Region 

4.  Solid-Liquid  Interface  Boundary  Layer 
Region 

5.  Isothermal  Inviscid  Core 


Figure  1.1  General  Features  and  Various  Regions  of  the  Flow 
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2.  TECHNICAL  BACKGROUND 
2.1  Materials  Processing  Processes 

It  has  been  recognized  by  several  researchers  [1-14]  that  convection 
plays  a  major  role  in  arc,  laser,  and  electron  beam  weld  pools.  Convection  is 
the  single  most  important  factor  influencing  the  geometry  of  the  pool 
including  pool  shape,  undercut,  and  ripples.  It  can  result  in  defects  such  as 
variable  penetration,  porosity,  and  lack  of  fusion.  It  is  also  primarily 
responsible  for  mixing,  thus,  affecting  the  composition  of  the  weld  pool. 

Fluid  flow  and  its  effect  on  weld  homogeneity  were  reported  by  Houldcroft 
[8]  and  Pumphrey  [9]  almost  three  decades  ago.  They  have  shown  that  for  a 
number  of  arc  weldings  with  a  filler  metal  of  different  composition  to  the 
base  metal  give  rise  to  substantially  uniform  weld  beads  which  imply  mixing  in 
the  weld  pool.  Also,  the  well-known  Shaeffler  diagram  [10]  implicitly  assumes 
complete  and  perfect  mixing  of  filler  metal  and  base  metal  plates  of  different 
composition,  and  is  found  to  work  out  well  in  practice. 

Apps  and  Milner  [11]  determined  the  "effective  thermal  conductivity"  of 
molten  lead  under  an  arc.  They  found  it  to  be  considerably  greater  than  that 
of  static  molten  metal,  showing  the  existence  of  convection  and  fluid  flow  in 
arc  weld  pools.  Both  Apps,  et  al.  [12]  and  Christensen,  et  al.  [13], 
attributed  the  shape  of  the  weld  pool,  at  least  partially,  to  the  motion  in 
the  molten  metal.  Bradstreet  |14j  reported  the  general  patterns  for  the  fluid 
flow  within  the  weld  pool  and  its  effect  on  bead  formation. 

Fluid  flow  in  weld  pools  also  affects  the  mixing  [15-19]  of  the  reacting 
products  and  the  solidifying  structure  [20,21].  Even  relatively  mild  con¬ 
vection,  as  experienced  in  ingot  casting,  affects  the  solidifying  structure 
[22].  Chande  and  Mazumder  [15],  observed  that  an  effective  diffusivity  of 
10®cm^/sec  (10  orders  of  magnitude  higher  than  the  actual  diffusivity)  is 
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required  to  explain  the  experimentally  obtained  uniform  solute  redistribution 
during  Laser  Surface  Alloying  (LSA)  on  the  basis  of  molecular  diffusion.  The 
implication  as  pointed  out  by  Chande  and  Mazumder  is  that  the  fluid  flow 
within  the  molten  pool  plays  a  significant  role  in  the  mechanism  of  solute 
redistribution. 

2.2  Surface  Tension  Gradient  Driven  Flow 

Surface  tension  gradient  driven  flow  has  been  identified  to  be  respon¬ 
sible  for  ripple  formation  in  the  weld  (23-26).  Copley,  et  al.  (23-25),  dis¬ 
cussed,  qualitatively,  the  mechanism  of  ripple  formations  by  surface  tension 
gradient  driven  flow.  Anthony  and  Cline  (26)  developed  a  steady  state  one¬ 
dimensional  model  for  the  fluid  flow  and  predicted  the  spacing  and  height  of 
the  surface  ripple.  The  momentum  equation  is  solved  based  on  the  temperature 
field  given  by  the  conduction  equation.  Therefore,  convective  heat  transfer 
within  the  molten  region  was  not  considered.  Oreper,  et  al.  (27),  did  a 
numerical  solution  for  the  convective  heat  transfer  within  the  molten  pool 
during  arc  welding.  Buoyancy,  electromagnetic,  and  surface  tension  forces 
were  considered.  It  was  found  that  the  surface  tension  gradient  is  the 
dominant  factor  in  many  cases.  The  model  has  the  inconsistency  that  the 
solid-liquid  interface,  which  is  not  known  a  priori,  was  being  specified. 
Heat  transfer  and  fluid  flow  during  electron  beam  welding,  when  keyholing 
phenomenon  occurs,  was  studied  by  W.  Geidt  (28).  In  such  a  case,  a  thin  layer 
of  molten  metal  is  found  between  the  plasma,  which  is  formed  by  the  keyholing 
effect,  and  the  solid  phase.  Within  the  framework  of  boundary  layer 
assumption,  he  obtained  a  solution  for  the  fluid  flow  and  the  heat  transfer. 
Surface  tension  gradient  driven  flow  also  arises  in  zone  melting  crystal 
growth.  A  general  formulation  and  dimensionless  governing  parameters  were 
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derived  by  Ostrach  [29].  An  extensive  review  of  literatures  in  this  area  can 


be  found  in  Ref.  [29). 


2.3  Multidimensional  Phase-Change  Problem 


The  problem  of  multidimensional  phase-change  heat  transfer  is  a  very  com¬ 


plicated  one.  It  involves  the  coupling  of  the  energy  equations  in  different 


phases  through  the  interface  which  in  most'cases  is  irregular.  Analytical  and 


semi-analytical  methods  have  a  very  limited  range  of  applicability  [30). 


Numerical  solutions  using  finite  difference  or  finite  elements  are  more 


versatile.  A  detailed  review  of  the  numerical  methods  has  been  published  by 


Shamsundar  [31].  He  compared  various  numerical  methods  of  solving  multi¬ 


dimensional  phase-change  problem.  Various  methods  of  solving  multidimensional 


phase-change  heat  transfer  were  reviewed  by  Viskanta  [30].  The  available 


methods  can  be  broadly  divided  into  two  groups,  (1)  temperature  as  the 


primitive  variable,  and  (2)  enthalpy  as  the  primitive  variable.  Large  number 


of  numerical  works  are  based  on  finite  difference  formulation,  although  finite 


element  seems  to  be  able  to  treat  the  irregular  boundaries  which  often  exist 


in  multidimensional  phase-change  problem.  A  general  finite  element  solution 


method  has  been  developed  by  Lynch  and  O'Neil  [32).  The  method  has  been 


applied  to  a  one-dimensional  phase-change  heat  transfer  problem  only. 


The  finite  difference  schemes  are  much  better  developed.  The  temperature 


based  method  can  be  further  divided  into  four  categories  [30]:  (1)  the 


explicit  finite  difference  scheme,  (2)  the  implicit  finite  difference  scheme. 


(3)  the  moving  boundary  immobiliration  method,  and  (4)  the  isotherm  migration 


method.  The  enthalpy  based  methods,  on  the  other  hand,  have  their  merits.  A 


detailed  discussion  of  the  advantages  and  disadvantages  of  the  enthalpy-based 


model  can  be  found  in  Ref.  (31).  Brief  descriptions  of  various  numerical 
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methods  (both  finite  difference  and  finite  element)  and  a  long  list  of 
references  can  be  found  in  the  review  by  Viskanta  [30]. 

Multidimensional  phase-change  with  convection  is  also  discussed  in 
[301.  In  addition  to  the  complication  of  coupling  of  the  energy  equations  of 
the  different  phases  through  the  interface,  there  is  the  coupling  of  energy, 
momentum,  and  continuity  equations  within  the  liquid  phase.  Some  works  have 
been  done  for  the  case  of  phase-change  with  natural  convection  in  the  liquid 
phase  [30]. 

In  order  to  achieve  a  quantitative  understanding  of  the  effects  of 
convection  on  pool  shape  and  mixing,  mathematical  modeling  of  the  flow  in  the 
pool  is  needed.  In  an  arc  weld  pool  both  surface  tension  gradient  driven  flow 
and  electromagnetic  field  driven  flow  play  a  role,  whereas  in  laser  welding 
only  the  flow  driven  by  surface  tension  variation  is  important.  Therefore, 
laser  melted  pool  is  an  ideal  candidate  to  study  the  affect  of  the  flow  driven 
by  surface  tension  variation  separated  from  electromagnetic  field  driven 
flow.  Recently,  it  has  also  been  recognized  that  at  moderate  current  (up  to 
150  ampere)  surface  tension  plays  a  major  role  even  in  arc  melted  pools 
[27).  In  view  of  the  foregoing,  convection  in  laser  melted  pools  is  analyzed. 

The  present  report  consists  of  three  models  of  the  thermocapillary 
convection  in  the  melt  pool.  The  first  model  is  a  perturbation  solution  of 
the  full  three-dimensional  problem.  As  it  turns  out  the  recirculatory  flow 
due  to  thermocapillarity  is  much  larger  than  the  scanning  speed.  Perturbation 
solution  is  therefore  sought.  This  approach  has  the  advantage  of  reducing  the 
original  three-dimensional  problem  to  two-dimensional.  However,  this  solution 
is  limited  to  lower  scanning  speed.  The  second  model  is  the  stagnation 
flow.  It  can  be  observed  from  the  results  of  the  perturbation  model  (or  the 
full  three-dimensional  model)  that  near  the  center  right  underneath  the 
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beam.  The  flow  has  a  great  similarity  to  stagnation  flow  except  the  motion  is 
now  driven  by  the  surface  shear  not  the  impinging  flow.  This  provides  a  clear 
picture  of  the  structures  of  the  fluid  flow  and  heat  transfer  in  this 
stagnation  region.  Quantities  such  as  maximum  temperature,  surface  velocity, 
viscous  and  thermal  boundary  layers  can  be  estimated  from  the  derived 
algebraic  equations.  The  final  model  is  a  full  three-dimensional  calculation 
of  the  Navier-Stokes  and  energy  equations.  This  model  also  has  the  capability 
to  treat  free  surface  flow. 

In  the  present  work,  we  have  limited  ourselves  to  the  problems  of  the 
fluid  flow  and  the  heat  transfer  in  the  molten  and  solid  regions.  The  laser 
beam  is  modeled  as  a  surface  heat  source  so  that  the  keyholing  phenomenon  is 
not  considered.  The  objective  of  this  research  is  to  understand  the  physics 
of  the  process--the  mechanism  and  the  coupling  of  the  fluid  flow  and  the  heat 
transfer.  The  effects  of  the  fluid  flow  and  the  heat  transfer  on  pool  char¬ 
acteristics,  such  as  pool  shape,  cooling  rate,  recirculating  flow,  and 
temperature,  are  quantified. 
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3.  FORMULATION 

3.1  Perturbation  Model 

3.1.1  Introduction 

In  this  section,  the  perturbation  model  of  the  three-dimensional 
problem  is  presented.  It  was  found  that  the  recirculating  velocity  is  much 
higher  than  the  scanning  velocity.  Therefore,  the  convection  due  to  the  re¬ 
circulating  velocity  is  much  larger  than  that  due  to  the  scanning  velocity. 
In  view  of  the  above,  a  perturbation  solution  is  therefore  attempted.  The 
basic  solution  corresponds  to  the  stationary  spot  source,  and  the  perturbation 
is  based  on  a  small  scanning  velocity.  The  advantage  of  seeking  a  pertur¬ 
bation  solution,  as  it  turns  out,  is  that  the  three-dimensional  problem  is 

solved  by  solving  two  sets  of  two-dimensional  equations.  The  formulation  is 
as  follows.  The  non-dimensional  form  is  first  derived.  The  velocity  field  is 
next  decomposed  into  two  parts,  the  first  part  is  the  recirculation  and  the 
second  is  the  scanning  velocity.  An  asymptotic  series  expansion  is  sought. 

The  zeroth  order  terms  correspond  to  the  axisymmetric  solution.  The  first 

order  terms  are  the  perturbations  due  to  low  scanning  speed.  The  angular 

dependence  of  the  perturbation  terms  are  separable  thereby  reducing  the  three- 
dimensional  problem  to  two-dimensional. 


3.1.2  Formulation 

The  final  dimensionless  governing  equations  are  given  below,  for 
detail  deriviation  refer  to  Ref.  [42]. 
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the  energy  equation  in  the  solid  region: 
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The  appropriate  boundary  conditions  are,  on  the  surface, 
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The  governing  equations  for  the  perturbation  are: 
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0  -momentum  equation: 
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z-momentum  equation: 
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the  energy  equation  within  the  molten  region: 
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the  energy  equation  in  the  solid  region: 
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the  boundary  conditions  become, 
on  the  free  surface, 
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along  the  solid-liquid  interface. 
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and  finally  at  a  distance  away. 
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RF=^(1^)^  (3.1.19) 

The  three-dimensional  problem  is  therefore  represented  by  two  sets  of  two- 
dimensional  equations.  The  basic  solution  is  described  by  the  system 
Eqs.  (3.1.2)  through  (3.1.9),  and  the  perturbation  solution  is  described  by 
the  system  Eqs.  (3.1.10)  through  (3.1.18).  However,  the  perturbation  model  is 
limited  to  slow  scanning  speed. 

Numerical  solutions  are  sought  for  the  two  sets  of  two-dimensional 
equations  of  the  perturbation  model.  Finite  difference  equations,  central 
difference  for  the  diffusion  terms  and  upwind  difference  for  the  convective 
terms,  governing  the  velocities  and  temperature  are  derived.  The  staggered 
grid  is  used.  The  resultant  non-linear  algebraic  equations  are  solved  by 
standard  Alternating  Direction  Iteration  (ADI)  method  [34]. 

The  computational  procedure  is  as  follows.  For  a  certain  temperature 
field  where  both  solid  and  liquid  region  exist,  the  interface  is  first  de¬ 
termined.  The  velocities  within  the  molten  pool  are  iterated  for  a  prescribed 
number  of  iterations.  Using  this  updated  velocity  field,  the  energy  equations 
are  iterated  next.  Thus,  the  temperature  field  for  the  next  global  iteration 
is  obtained.  The  initial  guesses  are  the  steady  state  conduction  temperature 
and  zero  velocity  for  the  basic  solution,  and  zero  temperature  and  velocity 
for  the  perturbation  solution.  The  convergent  solutions  are  obtained  when  the 
average  residual  of  each  equation  is  smaller  than  a  specified  small  number, 
typically  10"^. 

3.1.3  Results  and  Discussion 

Numerical  solutions  for  nine  cases  have  been  obtained.  The  values 
of  the  governing  parameters  are  tabulated  in  Table  3.1.1.  In  order  to 
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visualize  this  three-dimensional  fluid  flow  and  convection,  the  results  are 
presented  in  the  familiar  cartesian  coordinates,  as  depicted  in  Fig.  3.1.1, 
even  though  the  calculations  were  done  in  cylindrical  coordinates.  Using 
0.0001  for  Uj/Uj.,  a  solution  can  be  obtained.  The  values  of  the  governing 
parameters  are  Case  #1  tabulated  in  Table  3.1.1.  The  following  graphs  were 
plotted  based  on  this  result  unless  specified  otherwise. 

a.  Surface  Temperature:  A  three-dimensional  plot  of  the  surface  temper¬ 
ature  is  presented  in  Fig.  3.1.1  for  Case  #1.  The  horizontal  plane  is  the 
surface  plane  and  vertical  axis  is  the  dimensionless  temperature.  The 
temperature  distribution  is  very  much  expected.  The  temperature  attains  its 
maximum  at  the  center  underneath  the  beam  and  decreases  radially  outward.  It 
is  interesting  to  note  that  the  temperature  gradient,  the  driving  force,  has 
three  distinct  regimes.  Its  magnitude  increases  from  zero  at  the  center  in 
the  radial  direction  roughly  until  the  edge  of  the  laser  beam,  it  then  de¬ 
creases  and  finally  increases  again  as  the  flow  approaches  the  edge  of  the 
molten  pool.  The  first  stage,  being  so  close  to  the  center,  the  increase  in 
temperature  gradient  is  caused  by  surface  heating.  The  fluid  flow  driven  by 
the  surface  tension  gradient  becomes  important  in  the  second  stage  where  it 
tends  to  smooth  out  the  temperature.  Finally,  the  flow  approaches  the  edge  of 
the  molten  pool,  and  turns  around.  This  phenomenon  is  very  similar  to  the 
heat  transfer  near  the  stagnation  point,  and  therefore  gives  rise  to  high 
temperature  gradient. 

b.  Velocity  Field:  The  variation  of  surface  temperature  due  to  surface 
heating  induces  a  surface  tension  gradient  thus  pulling  the  molten  materials 
radially  outward.  The  velocity  field  on  the  surface  of  the  molten  pool  is 


plotted  in  Fig.  3.1.2.  The  velocity  field  in  the  vertical  plane  at  =  0  is 
presented  in  Fig.  3.1.3.  The  general  pattern  of  the  flow  field  is  that  the 
molten  materials  are  going  radially  outward  on  the  surface.  As  the  flow  ap¬ 
proaches  the  edge  of  the  molten  pool,  it  goes  down  and  turns  around.  It  then 
moves  to  the  center  and  comes  up  to  complete  the  recirculation  pattern.  In 
addition  to  this  recirculation,  there  is  also  the  motion  of  the  workpiece. 
Therefore,  material  enters  the  molten  pool  from  the  front  part  of  the  pool, 
goes  through  the  recirculation  pattern,  ultimately  resolidifying  on  the 
trailing  edge. 


c.  Solid  Redistribution:  In  order  to  help  the  reader  to  obtain  a  quali¬ 
tative  understanding  of  the  mechanism  of  solute  redistribution,  a  particle 
trajectory  is  plotted  in  Fig.  3.1.4.  This  trajectory  is  presented  in  three 
views;  front,  side,  and  top.  The  initial  position  is  (x^,  x^,  x|)  =  (-2.97, 
0.423,  0.01);  a  point  very  close  to  the  surface  and  before  the  molten  pool. 
The  recirculating  pattern  of  the  particle  can  be  clearly  observed.  This 
recirculating  pattern  implies  that  a  particle  will  travel  through  a  rather 
long  path  before  it  freezes  into  the  resolidifying  surface.  Consequently,  the 
molten  materials  can  be  well  mixed.  The  composition  within  the  melt  pool  is 
therefore  uniform. 


d.  Temperature  Field  within  the  Molten  Pool;  The  isotherms  within  the 
molten  pool  of  the  vertical  planes  at  x|  =  0  and  x|  =  0  are  plotted  in 
Figs.  3.1.5  and  3.1.6,  respectively.  Xj  is  the  scanning  direction.  The 
vertical  plant  at  x^  =  0  therefore  represents  the  cross-section  of  the  molten 
pool  directly  under  the  laser  beam.  It  can  be  observed  that  the  isotherms  in 
the  vicinity  of  the  solid-liquid  interface  are  compressed  near  the  surface  of 
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the  pool  and  stretched  apart  near  the  bottom  due  to  existence  of  the  flow 

S  field.  This  distortion  is  due  to  recirculation  of  molten  material.  The 

vertical  plane  =  0  is  the  longitudinal  section  of  the  molten  pool  along  the 
center  line.  The  asymmetry  of  the  isotherms  and  the  sol  id- liquid  interface 
can  be  observed.  This  distortion  is  due  to  the  motion  of  the  workpiece. 
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e.  Cooling  Rate:  To  visualize  the  heating  and  cooling  of  the  materials 
as  they  scan  pass  underneath  the  laser,  the  temperature  versus  the  scanning 
axis  are  plotted  in  Figs.  3.1.7  and  3.1.8.  In  Fig.  3.1.7,  the  family  of  the 
curves  represents  the  temperature  on  the  surface  plane  (x^  =  0).  Each  curve 
corresponds  to  a  different  value  of  x|.  On  the  other  hand,  the  curves  plotted 
in  Fig,  3.1.8  are  the  temperature  curves  in  the  vertical  plane  along  the 
centerline  (x|  =  0).  The  melting  temperature  is  also  plotted  in  both 
figures.  The  cooling  rate  right  after  resolidification  is  related  to  the 
microstructure.  The  higher  the  cooling  rate,  the  finer  is  the  micro¬ 
structure.  It  can  be  observed  from  the  previous  two  graphs  that  the  cooling 
rate  (i.e.  the  temperature  gradient)  decreases  from  maximum  the  centerline  to 
zero  at  the  edge  of  the  molten  pool  (Fig.  3.1.7)  and  also  decreases  from 
maximum  on  the  surface  to  minimum  at  the  bottom  of  the  molten  pool 
(Fig.  3.1.8).  Such  a  variation  in  cooling  causes  a  variation  of  the 
resolidified  microstructure.  Further  investigation  of  the  solidification 
mechanism  is  required  to  fully  predict  the  resulting  microstructure. 

f.  Pool  Shape;  The  molten  pool  is  greatly  altered  because  of  the 
presence  of  the  fluid  flow.  It  tends  to  make  the  molten  pool  wider  and 
shallower.  This  effect  is  shown  clearly  in  Fig.  3.1.9,  where  the  solid-liquid 
interfaces  at  the  vertical  section  (x|  =  0.0)  for  both  pure  conduction  and 
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convection  are  plotted.  The  convection  cases  are  tabulated  in  Table  3.1.1 
(Case  #1,4,  and  5).  The  Prandtl  number  is  the  only  changing  parameter  among 


these  cases.  An  increase  in  the  Prandtl  number  signifies  an  increased  in 
convection,  thus,  making  the  pool  wider  and  shallower. 

3.1.4  Conclusions 

A  perturbation  model  of  the  fluid  flow  and  heat  transfer  of  a  laser 
melted  pool  is  developed.  Because  of  the  fact  that  the  recirculating  flow  is 
much  greater  than  scanning  speed,  a  perturbation  is  sought.  The  basic 
solution  is  the  axisymmetric  case  and  the  perturbation  is  based  on  small 
scanning  speed.  The  advantage  is  that  the  three-dimensional  problem  is 
represented  by  two  sets  of  two-dimensional  governing  equations.  It  is  found 
that  the  presence  of  the  thermocapillary  convection,  the  physics  of  the 
process  changes  from  conduction  to  convection  dominated.  This  changes  the 
pool  geometry  dramatically,  resulting  in  up  to  a  150  percent  increase  in  the 
aspect  ratio  (width/depth)  as  compared  to  the  pure  conduction  case.  The 
temperature  gradient  on  the  resolidifying  surface,  which  is  directly  propor¬ 
tional  to  the  cooling  rate,  is  found  to  be  non-uniform.  This  in  term  gives 
arise  to  the  non-uniform  resolidified  microstructure.  The  mechanisms  of 
solute  redistribution  is  discussed  qualitatively.  It  is  found  that  a  particle 
recirculates  many,  many  times  in  the  molten  pool  before  it  resolidifies. 
Consequently,  the  solute  can  be  well  mixed. 
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Figure  3.1.2  Velocity  field  on  the  surface  plane  for 
Ma  =  1,500,  Pr  =  0.15,  T*  -  0.25  and 
RF  =  0.0  (perturbation  model) 
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igure  3.1.4  Trajectory  of  a  particle  as  it  enters,  recirculates 
and  ultimately  freezes.  The  trajectory  is  plotted  in 
front,  side,  and  top  view  in  (a),  (b),  and  (c), 
respectively.  Initial  position  is 
(X*  X*  X*)  =  (-2.97,  0423,  0.01)  Ma  =  1,500,  Pr  = 
0.15,  t*  =  0.25  and  RF  =  0.0  (perturbation  model) 


X^-AXIS 


Figure  3.1.7  Family  curves  represent  the  surface  temperature 
Each  curve  corresponds  to  a  different  value  of 


X*  for  Ma  =  1,500,  Pr  =  0.15,  T*  =  0.25, 
and  RF  *  0.0  (perturbation  modeT) 
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Figure  3.1.8  Family  curves  of  the  temperature  in  a -vertical  plane 
(X^  =  0).  Each  curve  corresponds  to  a  different 
value  of  X*  Ma  *  1,500,  Pr  =  0.15, 

T*  =  0.25and  RF  =  0.0  (perturbation  model) 
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Figure  3.1.9  Comparison  of  the  shape  of  the  molten  pool  for 
different  Prandtl  number  (perturbation  model) 
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Table  3.1.1  The  Numerical  Values  of  the  Dimensionless  Parameters 
of  the  Nine  Cases  of  the  Three-Dimensional  Model 


1 

2 

3 

4 

5 

Re 

10,000 

30,000 

30,000 

10,000 

10,000 

Pr 

0.15 

0.15 

0.015 

0.1 

0.2 

Ma 

1,500 

4,500 

450 

1,000 

2,000 

T* 

m 

0.25 

0.3 

0.3 

0.25 

0.25 

RF 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

7 

8 

9 

Re 

5,000 

5,000 

5,000 

10,000 

Pr 

0.10 

0.15 

0.2 

0.015 

Ma 

500 

750 

1,000 

150 

T* 

m 

0.25 

0.25 

0.25 

0.25 

RF 

0.0 

0.0 

0.0 

0.0 
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3.2  Stagnation  Flow  Analysis 
3.2.1  Introduction 

In  this  section,  a  generalized  solution  for  the  thermocapillary 
flow  and  convection  of  the  liquid  melt  under  an  intense,  non-uniform  heat 
source,  is  considered.  When  a  concentrated  heat  source  is  applied  to  a  free 
surface  of  a  molten  material,  the  resulting  temperature  distribution  causes  a 
non-uniform  surface  tension  distribution.  Being  a  decreasing  function  of 
temperature  for  most  liquid  metal,  the  surface  tension  gradient  will  therefore 
pull  the  liquid  from  the  warmer,  central  region  to  the  cooler,  outer  region  at 
the  surface.  This  shearing  motion  entrains  fluid  toward  the  surface  at  the 
center,  forming  viscous  and  thermal  boundary  layers  if  the  fluid  is  deep 
enough.  Such  flows  are  the  dominant  features  in  the  central  region  of  the 
melt  pools  during  laser  materials  processing  processes  such  as  laser  surface 
alloying  and  rapid-solidification  by  laser  melt  quenching  [15,27].  Similar 
flows,  with  a  different  geometry  are  also  found  in  crystal  growth  processes 
using  zone  melting  and  similar  methods  (29).  The  physical  model  is  shown  in 
Fig.  3.2.1.  Both  two-dimensional  and  axi symmetric  cases  are  considered.  The 
two-dimensional  case  represents  the  limiting  case  of  the  flows  encountered  in 
zone  melting  crystal  growth,  when  the  heat  source  dimension  is  small  with 
respect  to  the  radius.  The  axi symmetric  case  represents  the  limiting  case  of 
the  flows  encountered  in  welding,  laser  surface  alloying  and  rapid- 
solidification  by  laser  melt  quenching  when  the  scanning  velocity  is  small 
relative  to  the  velocity  of  the  thermocapillary  motion.  This  investigation  is 
to  focus  on  the  basic  features  of  the  flow  in  the  central  region  where  the 
structure  is  very  much  like  the  stagnation  flow.  The  theory  approximates  the 
surface  heat  flux  distribution  near  the  center  by  a  parabola,  and  invokes  the 
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boundary  layer  approximation  for  the  energy  equation,  but  not  for  the  momentum 
equation.  The  governing  coupled  partial  differential  equations  are  trans¬ 
formed  to  a  set  of  decoupled  ordinary  differential  equations.  Analytical 
solutions  are  found  for  some  conditions.  Numerical  solutions  are  sought  for 
other  conditions.  These  results  provide  the  basic  scaling  laws  for  thermocap¬ 
illary  convection  in  a  deep  melt  pool.  Explicit  formulas  are  obtained  for  the 
scalings  of  the  physical  quantities  such  as  velocity,  temperature,  surface 
temperature  gradient,  and  viscous  and  thermal  boundary  layers.  Because  of  the 
intrinsic  coupling  of  the  momentum  and  energy  equations  of  thermocapillary 
convection,  these  scalings  depend  on  the  Prandtl  number.  This  dependence  is 
obtained  by  computations  of  the  energy  equations  for  various  Prandtl  numbers 
near  unity.  The  limiting  cases  of  large  and  small  Prandtl  numbers  are  also 
analyzed.  The  results  complete  the  understanding  of  scaling  laws  dependence 
on  the  Prandtl  numbers.  These  asymptotic  analyses  also  yield  explicit  and 
concise  solutions  and  scalings  within  the  respective  ranges.  It  is  found  that 
the  asymptotic  results  give  good  approximation  for  moderate  values  of  the 
Prandtl  numbers  as  compared  to  the  numerical  results  of  the  case  Pr  -  0(1). 

3.2.2  Formulation 

The  governing  equations  are  given  below,  for  detail  deriviation 
refer  to  Ref.  [42]. 

Pr  -  0(1) 

fj"  f  (m-i)  r;  -  =  0 


(3.2.1) 
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g"  +  Pr  F,  g;  =  0 


gj  +  Pr  (Fj  g|  -  2  F|  g|  =  0 


the  boundary  conditions  are. 


0;  F^  =  0,  F''  =  -1,  g;  =  g^  =  -1 


z  -  ®  F^  -  §0  =  g^  =  0. 


Pr  -  0 


g  °  +  F  g  °  =  0, 


(3.2.2) 


(3.2.3) 


(3.2.4) 


(3.2.5) 


(3.2.6) 


where  6(z°)  is  the  Dirac  delta  function  and  =  f|^(”)  is  a  constant 
listed  in  Table  3.2.1. 


the  boundary  conditions  are. 


2°  =  0,  g^,°'  =  gj°'  =  -1, 


(3.2.7) 


•  o  ^Tsnx.fi^.ro  o  o  4-*V^  ^  . 
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where  =  Fj^(O)  is  a  constant  tabulated  In  Table  3.2.1. 


the  boundary  conditions  are, 

2^  =0;  §0  =  gj  =  -1. 

♦  «.;  §0  =  g|  *  0.  (3.2.10) 
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3.2.3  Results  and  Discussion 

a.  Pr  ~  0  (1):  The  universal  stream  function  f^,  velocity  pro¬ 
file  f^,  and  shear  stress  F'^  for  both  two-dimensional  and  axisymmetric  cases 
are  plotted  in  Fig.  3.2.2.  The  behavior  of  the  stream  function  is  quite  dif¬ 
ferent  from  the  conventional  stagnation  flow.  In  the  case  of  thermocapillary 
flow,  the  fluid  flow  is  driven  by  the  surface  tension  gradient  on  the  free 
surface.  Consequently,  the  shear  stress  is  negative  with  maximum  magnitude  on 
the  free  surface  and  diffuses  into  the  liquid.  The  velocity,  u^,  is  non-zero 
on  the  surface.  In  fact,  it  is  maximum  here  and  decreases  exponentially  with 
the  depth.  The  scaled  temperature  functions  g^  and  g^  for  various  Prandtl 
numbers  are  plotted  in  Figs.  3.2.3  through  3.2.10  for  two-dimensional  and 
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ax i symmetric  cases.  The  scaled  temperature  functions  decrease  very  rapidly 
from  their  maximum  values  on  the  surface.  The  maximum  values  which  exist  on 
the  surface  increase  as  the  Prandtl  number  decreases.  This  is  because  the 
smaller  the  Prandtl  number,  the  smaller  is  the  convection,  less  heat  being 
brought  sideways,  hence  the  higher  is  the  temperature. 

The  decoupling  of  the  system  by  simple  rescaling  means  that  the  coupling 
effect  of  the  heat  transfer  and  momentum  is  by  a  multiplication  factor.  This 
can  be  seen  easily  by  writing  the  viscous  boundary  layer  scale,  velocity,  the 
maximum  temperature,  and  the  heat  transfer  coefficient  explicitly  as: 


\  ^lo 


^lo 


'  ‘'l'’  5lo'’  >'*  fl<^> 


_  0  /  kuv  z  Z  /n\ 

’  r  <7^57)  ’10  9o(°> 


q  g  (0)  ^ 

^0^0'  ,  uv 

qfgio^  Jx't' 


(3.2.11) 


(3.2.12) 


(3.2.13) 


g^^  appears  in  all  these  quantities,  g^^  is  a  measure  of  the  surface  temper¬ 
ature  non-uniformity.  In  particular,  it  is  the  second  derivative  (with 
respect  to  r)  of  the  surface  temperature.  This  is  introduced  into  the  non- 
dimentionalization  to  decouple  the  system.  Thus  the  coupling  of  the  energy 
and  momentum  equations  is  solely  embodied  in  this  term,  g^^^  is  a  function  of 
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the  Prandtl  number  and  is  graphically  presented  in  Fig.  3.2.11  for  both  two¬ 


dimensional  and  axisymmetric  cases.  gQ(0)  denotes  the  dimensionless  surface 


temperature  at  the  center  (r  =  0)  and  therefore  is  directly  related  to  the 


maximum  temperature.  Naturally,  gQ(0)is  also  a  function  of  the  Prandtl  number 


and  is  plotted  in  Fig.  3.2.12.  The  surface  temperature  gradient  can  be 


obtained  by  taking  the  derivative  of  the  temperature  with  respect  to  r  and 


then  by  substituting  in  the  scalings.  It  is  found  to  be: 


ill  =  q  3/4  ^ 

I-  ain  ^ 


k  ^2Y'q^^  ^lO 


(3.2.14) 


i  iO  2k^Y' 


The  dependence  of  all  the  scalings  on  the  temperature  coefficient  of  surface 


tension,  and  other  thermophysical  properties  are  expressed  explicitly  in 


Eqs.  (3.2.1)  through  (3.2.4).  AS  mentioned  earlier,  g.  appears  in  all  of  the 


scalings,  and  it  is  a  function  of  the  Prandtl  number.  The  complete 


understanding  of  the  scalings  required  the  asymptotic  analyses  of  small  and 


large  Prandtl  numbers  which  are  presented  in  the  following  two  sections, 


b.  Pr  -►  0:  The  dimensionless  functions  g°(z°)  and  g°(z°)  are  plotted  in 


Fig.  3.2.13.  This  limiting  analysis  provides  explicit  formulas  for  the 


Prandtl  number  effect  on  the  scaling  laws.  The  scaling  introduced  for  gj^(z) 


[421  implies  that  the  viscous  length  scale  is: 


(3.2.15) 


..'i 


>  •  -  ’ 


This  viscous  length  sdale  can  be  rewritten  as; 


( 


2r  3 


2y' 


-) 


(3.2.16) 


3/4 

Thus,  the  viscous  boundary  layer  is  proportional  to  v  .  The  ratio  of 
the  thermal  length  scale  and  viscous  length  scale  is. 


6  /6  =  Pr 

<  V 


(3.2.17) 


so  that 


7  _  /  kyv  ,1/4  -o  ^  p^-3/4 

^  -  ^2Y'qi^  ^lo 


2  CJ  y'v 


(3.2.18) 


The  same  Prandtl  number  dependence  is  found  for  the  temperature.  In 
particular,  the  maximum  temperature  difference  which  occurs  on  the  surface  at 
the  center  can  be  obtained.  Substituting  the  scaling  for  g^,  we  have. 


at 


max 


^0  /  kpv  xl/4 
k  ^2Y'qi' 


gO  (0) 
g°(0)^/^ 


(3.2.19) 


The  driving  force  of  the  flow  field,  the  surface  temperature  can  be  shown  to 
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-  ,^-3/4  .  , 

,  „  3/4  =0  ,  1  ,'• 

=  2  ■’l  9lo  <2  ,  ,3  > 

2  p  Y  Cp  V 


(3.2.20) 


The  velocity  in  the  radial  direction  is 


—  r  fl(2) 
1/0  ^ 


50  -1/^  1/2 

^2Y'q^^  ^lo 


=  "l'*  9?o'’  9?o'*  ■  ^ 


3.2.21) 


The  vertical  velocity,  u,,  can  be  derived  from  the  continuity  equation: 


-  (fn-U  V 

,  kyv  .1/4  -0 
2y  qj^  lo 


:i74^<^> 


(3.2.22) 


<"-1)  "i'*  9?o'’  fi(^> 

P  P 


The  shear  stress  can  be  obtained  by  differentiating  u^  with  respect  to  z  in 
Eq.  (3.2.11) 


(3.2.23) 


=  q 


3/4 

1 


3/4 


r 


which  appears  in  the  above  equations,  is  a  constant  tabulated  in 
Table  3.2.2  for  both  two-dimensional  and  axisymmetric  cases. 


c.  Pr  -  The  dimensionless  functions  9^(2^)  and  gj(z^)  for  the  two- 
dimensional  and  axisymmetric  cases  are  plotted  in  Fig.  3.2.14.  The  solutions 
behave  as  expected,  attending  their  maximum  on  the  surface  and  decreasing  very 
rapidly  downward. 

Similarly,  this  asymptotic  limit  provides  the  explicit  formula  for  the 
effect  of  the  Prandtl  number  on  the  scaling  laws  for  large  Pr.  The  viscous 
length  scale  is; 


kuv  \ 
2Y'qi^ 


1/4 


Prl/8 


(3.2.24) 


=  a  qi 
'’l 


-H  k  p  C  1/8 

( - 


The  ratio  of  the  thermal  length  scale  and  viscous  length  scale  is; 
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(3.2.25) 


so  that 


kyiNJ  \ 
2Y'qi^ 


p,-3/8 
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f'-l 
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1/8 
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(3.2.26) 


It  is  interesting  to  compare  this  high  Prandtl  number  limit  of 
the  s  /6  ratio  of  the  present  model  to  that  of  the  common  forced  con- 

<  V 

vection.  The  limit  for  common  forced  convection  case  is  6  /6  <==  Pr'^/^.  This 

K  \) 

difference  is  a  result  of  the  non-zero  velocity  on  the  surface  so  that  the 
stream  function  behaves  linearly  in  z  for  small  z.  On  the  other  hand,  the 
common  forced  convection  problem  has  the  no-slip  condition  so  that  the  stream 
function  behaves  quadratical ly  in  z  for  small  z.  The  maximum  temperature 
difference  is  found  to  be, 
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is  a  constant  tabulated  in  Table  3.3.2  for  both  two-dimensional  and  axi- 
symmetric  cases,  for  the  whole  range  of  Pr,  which  is  necessary  for  the 
scaling  laws,  is  also  presented  graphically  in  Fig.  3.2.15  for  both  two- 
dimensional  and  axisymmetric  cases.  The  maximum  temperature  as  a  function  of 
the  Prandtl  number,  including  both  high  and  low  Pr  limits,  is  plotted  in 
Fig.  3.2.16.  The  asymptotic  limits  are  reflected  in  the  graph  (log- log  plot) 
as  two  straight  lines  of  slopes,  -1/2  and  -1,  respectively.  It  can  be  ob¬ 
served  that  the  two  asymptotic  limits  give  with  reasonable  predictions  of 
these  quantities.  For  Pr  <  -.15,  the  maximum  errors  of  the  low  Pr  limit 
prediction  for  both  9^(0)  and  §^^(0)  are  less  than  10  percent.  For  Pr  >  10.0, 
the  maximum  errors  of  the  high  Pr  limit  prediction  for  both  9^(0)  and  9^(0) 
are  also  less  than  10  percent. 

d.  Comparison  to  Numerical  Solution:  The  sta9nation  flow  analysis 
presented  in  this  section  is  the  dominant  feature  of  the  thermocapillary 
convection  in  melt  pool.  This  model  represents  a  common  feature  of  thermo¬ 
capillary  convection  in  the  melt  pool  in  the  vicinity  of  the  central  re9ion 
underneath  the  heat  source.  In  order  to  examine  the  validity  of  the  model, 
this  stagnation  flow  analysis  results  are  compared  to  the  numerical  results  of 
the  axisymmetric  case  presented  in  the  last  chapter. 

i.  Velocity  Profile 

The  universal  velocity  profile  of  the  stagnation  flow  analysis  is 
plotted  in  Fig.  3.2.17.  The  velocity  profile  obtained  from  the  perturbation 
model  (case  #2  in  Table  3.1.1)  are  plotted  also  in  Fig.  3.2.17  to  provide  the 
comparisons  among  themselves  as  well  as  to  the  stagnation  flow  analysis 
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result.  The  comparison  shows  that  the  surface  velocities  predicted  by  the  two 
theories  are  within  10  percent.  These  predicted  velocity  profiles  are  also 
similar.  The  returned  flow  velocity  at  the  bottom  part  of  the  pool  which  is 
neglected  in  the  stagnation  flow  analysis  is  predicted  to  be  of  one  order  of 
magnitude  less  than  the  surface  velocity  by  the  numerical  solution.  The 
stagnation  flow  analysis  predicted  that  the  velocity  profile  is  directly 
proportional  to  r.  It  can  be  observed  that  the  curves  f'"'  versus  z  for 
different  values  of  r*  of  the  numerical  solutions  are  essentially  one  curve. 

ii.  Centerline  Temperature 

The  centerline  temperature  predicted  by  the  stagnation  flow  analysis 
(i.e.  gQ(z)  versus  2  for  Pr  =  0.15)  is  plotted  in  Fig.  3.2.18.  This  center 
line  temperature  of  case  #2  (Table  3.1.1)  is  also  plotted  in  Fig.  3.2.18  for 
comparison  of  the  two  predictions.  The  stagnation  flow  analysis  predicts  a 
high  surface  temperature,  10  percent  higher  than  the  numerical  solution.  This 
is  expected  since  the  boundary  layer  approximation  is  invoked  for  the  energy 
equation  so  that  the  radial  conduction  is  not  taken  into  consideration.  The 
center  line  temperature  drops  more  rapidly  in  the  case  of  the  stagnation  flow 
analysis.  This  is  because  of  the  fact  that  the  stagnation  flow  analysis  is 
based  on  the  boundary  layer  approximation  of  the  energy  equation  and  the 
assumption  of  semi-infinite  domain.  The  thermocapillary  convection,  in  the 
case  of  stagnation  flow  analysis,  since  there  are  no  return  flow,  carries  all 
the  heat  from  the  heat  source  sideways.  However,  in  the  case  of  the  numerical 
solution,  there  is  a  return  flow  which  convects  part  of  the  heat  back. 


iii.  Maximum  Temperature 

The  maximum  temperature,  which  occurs  at  r*  =  z*  =  0,  of  all  the 
numerical  solution  cases  and  stagnation  flow  analysis  are  plotted  in 
Fig.  3.2.19.  The  stagnation  flow  analysis  is  predicting  the  maximum  temper¬ 
ature  quite  accurately  (within  10%)  for  Pr  =  0.10,  0.15,  and  0.2  However,  the 
stagnation  flow  analysis  is  predicting  a  much  higher  temperature  for  small 
Prandtl  number  (Pr  =  0.015).  This  is  due  to  the  boundary  layer  approximation 
on  the  energy  equation.  In  essence,  the  radial  conduction  is  neglected. 

3.2.4  Conclusions 

A  generalized  solution  for  the  thermocapillary  convection  in  the 
central  region  of  the  liquid  melt  pool  directly  underneath  an  intense  non- 
uniform  heat  source  is  obtained  for  both  two-dimensional  and  axisymmetric 
cases.  The  system  of  coupled  partial  differential  equations  is  transformed 
into  a  set  of  decoupled  ordinary  differential  equations.  By  doing  so,  a 
universal  stream  function  (velocity  profile)  is  obtained.  The  rescaled 
temperature  field  is  solved  for  various  Prandtl  numbers.  The  effects’  of  the 
coupling  of  the  temperature  and  the  velocity  fields  are  characterized  by 
scalings.  The  term  g^^,  appeared  in  all  of  these  scalings,  represents  the 
leading  order  term  of  the  temperature  non-uniformity  on  the  surface  and  varies 
with  respect  to  the  Prandtl  numbers.  Quantitative  dependence  of  g^^  on  the 
Prandtl  numbers  is  obtained  and  it  is  presented  graphically.  This  analysis 
also  yields  explicit  formulas  of  physical  quantities  such  as  maximum  temper¬ 
ature,  velocity,  shear  stress,  and  viscous  boundary  layer  length  scale.  These 
expressions  are  presented  in  Eqs.  (3.2.1)  through  (3.2.4).  Their  dependences 
on  thermophysical  properties  are  explicitly  defined.  The  two  asymptotic 


limits  of  small  and  large  Prandtl  numbers  provide  concise  formulas  for  these 
dependences  in  the  respective  ranges.  These  two  asymptotic  limits  also  yield 
explicit,  simple,  and  concise  formulas  for  the  physical  quantities  such  as 
maximum  temperature,  second  derivative  of  the  surface  temperature,  velocity, 
shear  stress,  thermal  boundary  layer  length  scale,  and  viscous  boundary  layer 
length  scale.  These  expressions  are  presented  in  Eqs.  (3.2.5)  through 
(3.2.13)  and  Eqs.  (3.2.14)  through  (3.2.21)  for  small  and  large  Prandtl 
numbers,  respectively. 

The  results  of  the  present  stagnation  flow  analysis  are  also  compared 
with  the  numerical  solutions  of  the  base  flow  of  the  perturbation  model.  The 
results  of  the  numerical  solutions  are  plotted  in  the  stagnation  flow  analysis 
variables.  It  is  found  that  the  stagnation  flow  analysis  predictions  are 
close  to  the  numerical  solutions.  Thus,  the  simple  algebraic  formulas  of  the 
stagnation  flow  analysis  can  be  used  to  predict  various  physical  quantities 
such  as  the  maximum  temperature  and  the  velocity. 
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Graphs  of  g  versus  the  Depth  of  the  Two 
Dimensiona r Case  for  Various  Prandtl 
Numbers  (>  1) 


Graphs  of  g  versus  the  Depth  of  the  Two 
Oimensional°Case  for  Various  Prandtl 
Numbers  (<  1) 


Figure  3.2.8  Graphs  of  g  versus  the  Depth  of 
Case  for  Various  Prandtl  Numbers 
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Figure  3.2.15  Second  Coefficient  (Second  Derivative)  of  the 
Surface  Temperature  versus  the  Prandtl  Numbers 
(Stagnation  Flow  Analysis) 


Figure  3.2.16  Maximum  Temperature  versus  the  Prandti  Number 
(Stagnation  Flow  Analysis) 
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Figure  3.2.17  Comparison  of  the  Universal  Velocity  Profile  Fj 

for  the  Stagnation  Flow  Analysis  and  the  Numerical 

Solution  Results  Ma  =  1,500,  Pr  =  0.15, 

T*  =  0.25,  and  RF  =  0.0 
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Figure  3.2.19  Comparison  of  the  Maximum  Temperature  of  the 
Stagnation  Flow  Analysis  and  the  Numerical 
Solution  Results.  The  Numerical  Solution  Results 
are  Numbered  According  to  Table  5.1 
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Table  3.2.1  Asymptotic  Limits  of  the  Stream  Function  at  Small  and  Large  z 
Functional  Form  Plane  2-D  Axisymmetic 


z  ->■  0 
2  -►  00 


f  =  f  z 
1  lo 

i  “ 


f;  =  1.0 
lo 

f  =  1.0 
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f  =  0.7124 


Table  3.2.2  Asymptotic  Limits  of  the  Surface  Temperature  Coefficients 


Pr  -  0 

Plane 

2-0 
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Ax i symmetric 
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gj(0)  =  1.3946 
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3.3  Three  Dimensional  Numerical  Model 

3.3.1  Introduction 

The  majority  of  the  previous  work  dealing  with  surface  tension  gra¬ 
dient  driven  flows  was  based  on  one-dimensional  and  two-dimensional  assump¬ 
tions.  Because  of  the  complex  nature  of  the  melt  pool  phenomenon,  important 
physical  characteristics  may  be  unresolved  by  these  oversimplified  models. 

The  perturbation  solution  described  earlier  in  this  report  offers  an 
interim  solution  in  that  two  two-dimensional  solutions  are  used  to  construct 
three-dimensional  solutions  of  modest  scanning  velocities.  This  chapter 
describes  full  three-dimensional  numerical  solutions.  Our  program  is  a  three- 
dimensional  finite  difference  scheme  which  iteratively  evaluates  the  energy, 
momentum  and  continuity  (pressure  corrector)  equations  for  transient/  steady- 
state  cases.  For  programming  simplicity  as  well  as  easy  user  understanding,  a 
rectangular  coordinate  system  was  defined  which  uses  in-line  grid  and 
primitive  variable  formulations.  Assumptions  made  to  simplify  the  problem 
are:  1)  the  material  is  homogeneous  throughout,  2)  all  physical  properties, 
with  the  exception  of  viscosity  and  surface  tension,  remain  constant 
throughout  both  solid  and  liquid  phases,  3)  the  material  scanning  speed 
remains  constant,  4)  the  absorbed  surface  energy  flux  remains  constant,  5)  the 
top  surface  was  initially  flat  and,  6)  buoyancy,  evaporation  and  latent  heat 
effects  were  neglected. 

3.3.2  Formulation 

The  basic  governing  equations  in  indicial  notation  are. 


the  continuity  equation: 


kinematic  free  surface  condition. 
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dynamic  free  surface  condition, 
normal  stress  balance. 
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tangential  stress  balance. 
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heat  flux  condition. 


{ 


-q"  +  e  o  (T^  -  T^),  x^^  +  x/  <  r^^ 
4  4% 

e  o  (T  -  T,),  otherwise 

a 


(3.3.8) 


The  boundary  conditions  along  the  solid-liquid  interface  which  is  defined  by 
the  surface  at  melting  temperature. 


are  the  continuity  of  mass. 
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It  is  important  to  point  out  that  the  surface  which  describes  the  solid-liquid 
interface  is  not  known  a  priori.  It  is  part  of  the  problem  to  be  solved. 

The  formulation  is  completed  by  specifying  the  boundary  condition  far 
away  from  the  heat  source  to  be  unaffected  by  the  heat  source. 


hj,  Ix^l,  X3  *  »; 

T  -  (3.3.12) 

The  calculation  domain  is  illustrated  in  Figure  3.3.1.  A  right-handed 
rectangular  coordinate  system  is  imposed  on  an  elongated  slab,  the  origin  of 
the  coordinate  system  being  located  at  the  slab's  lower,  left,  front  corner. 
The  concentrated  heat  source  strikes  the  workpiece  on  the  top  surface,  the 
center  of  the  source  being  located  at  the  slab  mid-point  in  the  Y  direction 
and  approximately  one-third  of  the  slab  length  in  the  X  direction.  The  viewer 
and  the  heat  source  are  located  in  a  stationary  reference  frame  while  the 
workpiece  scans  into  and  out  of  the  "control  volume"  calculation  domain. 
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Because  the  material  is  assumed  to  be  scanning  in  the  positive  X  direction 
only,  symmetry  can  be  assumed  about  the  slab  mid-point  in  the  Y  direction. 
Due  to  this  symmetry  consideration,  the  actual  calculation  domain  is  reduced 
and  is  defined  by  the  bold,  solid  rectangular  box. 

Figure  3.2.2  shows  an  X-Z  cross-section  of  the  calculation  domain  located 
at  the  symmetry  plane  in  Y.  The  concentrated  heat  source  is  assumed  to  have  a 
Gaussian  energy  distribution  and  is,  again,  applied  at  the  top  surface  of  the 
calculation  domain.  Phantom  nodes  were  introduced  to  facilitate  the  impo¬ 
sition  of  desired  boundary  conditions.  A  Y-Z  cross-section  of  the  calculation 
domain  located  at  X  =  X  (source  center)  is  shown  in  Fig.  3.3.3.  As  seen,  the 
Gaussian  energy  distribution  is  centered  on  the  symmetry  plane;  the  actual 
calculation  domain  being  represented  by  the  bold,  solid  box. 

The  details  of  the  numerical  scheme  are  described  in  the  MS  Thesis  by 
R.  Zehr  (46).  It  is  essentially  a  point-iterative  computation  based  on  the 
finite  difference  formulation  of  the  governing  equations.  In  brief,  at  each 
iteration  step,  the  residuals,  representing  errors,  are  computed  for  the  three 
momentum  equations,  the  continuity  equation,  and  the  energy  equation.  These 
residuals  serve  as  indications  for  convergence,  as  well  as  the  bases  for 
computing  convections  for  each  of  the  velocities,  pressure,  and  temperature. 
The  convections  for  velocities  and  temperature  are  based  on  a  mixed  Gauss- 
Seidel  and  Jacobi  schemes,  proceeding  from  line-to-line  according  to  the 
Gauss-Seidel  scheme,  updating  values  as  each  line  is  completed.  With  each 
line,  however,  the  point-by-point  iteration  is  based  on  the  Jacobi  scheme 
without  updating  the  values.  This  is  doen  to  take  advantage  of  the 
vectorizability  of  the  Jacobi  scheme.  Pressure  corrections  are  based  on  the 
artificial  compressibility  concept  of  Chorin  [43],  using  a  form  for  implicit 
problems  derived  by  Mech  [44]. 
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JC*  For  the  calculations  involving  free  surface  deformation,  the  pressure 

|l  distribution  from  flow  calculations  and  that  from  the  surface  curvature  are 

compared  and  iterated  until  they  match.  For  these  and  other  details,  consult 
the  thesis  by  Zehr  [46]. 
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3.3.3  Numerical  Results 

a.  Introduction:  Numerical  solutions  were  obtained  for  six  test 
cases  in  which  steel  was  used  as  the  scanning  substrate.  The  governing  param¬ 
eters  are  listed  in  Table  3.3.1.  Table  3.3.2  shows  the  various  test  case 
powers  and  scanning  speeds  which  range  from  200  W  to  400  W,  and  0.01  m/s  to 
0.024  m/s,  respectively.  Selected  numerical  information  for  test  cases  I,  II, 
and  V  are  graphically  presented,  corresponding  to  Power  =  200  W,  =  0.010 
m/s;  Power  =  240  W,  =  0.012  m/s;  and  Power  =  400  W,  =  0.02  m/s, 
respectively.  It  should  be  noted  that  the  beam  center  (at  the  domain  surface) 
is  defined  as  (X  =  0,  Y  =  0,  Z  =  0)  for  graphing  purposes. 

b.  Temperature  Distributions:  Figures  3.3.4  through  3.3.12  illustrate 
temperature  isotherms  (contours)  from  various  cross-sections  (viewpoints)  of 
the  domain.  The  three  cross-sectional  views  are,  1)  X-Z  plane  with  Y  located 
at  the  symmetry  plane,  2)  Y-Z  plane  with  X  located  at  the  beam  center,  and  3) 
X-Y  plane  with  Z  located  at  the  material  surface.  As  seen,  the  temperature 
isotherms  compress  near  the  top  solid/liquid  interface  due  to  the  rapid 
impingement  of  fluid  (convecting  heat)  against  the  solid  substrate.  These 
isotherms  are  compressed  slightly  more  at  the  leading  edge  of  the  melt  pool  as 
compared  to  the  trailing  edge  due  to  the  scanning  velocity  of  the  substrate. 
The  gap  separating  the  isotherms  widens  near  the  bottom  of  the  melt  pool 
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indicating  slower  (convective)  velocities  in  this  region.  Figures  3.3.6  and 
3.3.10  illustrate  the  circular  nature  of  the  isotherms  near  the  (beam)  center 
which  transform  into  ellipses  due  to  the  influence  of  the  substrate  scanning 
speed  near  the  solid/liquid  interface.  This  axisymmetric  temperature  distri¬ 
bution  near  the  pool  center  supports  the  perturbation  approach  taken  by  Chan. 
[42]  Figures  3.3.7  and  3.3.11  represent  three-dimensional  surface  temperature- 
contours  which  vividly  illustrate  the  severe  temperature  gradients  beneath  the 
beam  and  at  the  solid/liquid  interface.  An  interesting  feature  is  the  dif¬ 
ference  between  the  upstream  and  downstream  dimensions  of  th  pool.  The  rela¬ 
tive  downstream  velocities  are  greater  than  the  relative  upstream  velocities 
causing  the  melt  pool  front  and  sides  (-X  and  ±Y,  respectively)  to  flatten. 
More  heat  is  convected  downstream  which  lengthens  (and  deepens)  the  melt  pool 
trailing  edge.  (See  Figures  3.3.4,  3.3.6,  3.3.8,  3.3.9,  3.3.10,  and  3.3.12.) 

c.  Velocity  Fields:  The  recirculating  velocity  fields  for  the  respective 
cross-sections  are  shown  in  Figures  3.3.13  through  3.3.17.  Fluid  is  driven 
radially  outward  along  the  surface  from  the  beam  center  until  it  reaches  the 
solid/liquid  interface.  At  this  point  the  flow  is  turned  downward,  circu¬ 
lating  along  the  pool  bottom  until  reaching  the  center  where  it  is  forced  up¬ 
wards  toward  the  surface.  This  recirculating  pattern  resembles  a  toroidal 
vortex.  The  large  radial  surface  velocities  are  a  direct  consequence  of  the 
dramatic  surface  temperature  gradients  shown  in  Figures  3.3.7  and  3.3.11. 
Numerical  results  indicate  that  thermocapillary  (surface)  velocities  are 
approximately  two  orders  of  magnitude  larger  than  the  substrate  scanning 
velocity  which  illustrates  the  importance  of  surface  tension  effects. 
Governing  dimensionless  parameters  and  corresponding  boundary  layer  thickness 
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estimates  are  listed  in  Table  3.3.3.  Reynolds  and  Peclet  numbers  are  based  on 
the  calculated  pool  length  (in  X)  and  the  maximum  surface  velocity  (U)  for 
each  test  case.  As  the  test  case  power  and  scanning  speed  increased  the  pool 
length  and  maximum  surface  velocity  also  increased.  This  increase  caused  the 
Reynolds  and  Peclet  numbers  to  increase  which  caused  the  viscous  and  thermal 
boundary  layers  to  grow  thinner. 

d.  Cooling  Rates:  To  illustrate  heating  and  cooling  of  the  substrate  as 
it  passes  beneath  the  laser,  temperature  versus  scanning  distance(X)  is 
plotted.  Figures  3.3.19  and  3.3.20  show  surface  temperature  versus 
distance(X)  for  various  locations(Y)  from  the  symmetry  plane  while  Figures 
3.3.19  and  3.3.21  illustrate  symmetry  plane  temperatures  versus  distance(X) 
for  various  locations(Z)  from  the  domain  surface.  The  substrate  is  observed 
to  increase  in  temperature  as  it  approaches  the  laser  and  decrease  in  temper¬ 
ature  after  passing  beneath  the  laser.  The  resolidified  substrate  micro¬ 
structure  is  controlled  by  the  magnitude  of  the  cooling  rate  at  any  given 
location;  the  greater  the  cooling  rate  the  finer  the  microstructure.  (The 
cooling  rate  is  directly  proportional  to  the  temperature  gradient.)  As 
illustrated,  the  most  severe  temperature  gradients  are  located  at  the 
(symmetry  plane)  surface  solid/liquid  interface,  becoming  less  severe  at  the 
melt  pool  side  edges  and  bottom.  These  results  suggest  that  the  resolidified 
substrate  microstructure  is  finest  at  the  symmetry  plane  surface  and  coarser 
at  the  side/bottom  edges.  This  is  contradictory  to  experimental  observation, 
however,  where  the  finest  microstructure  is  found  at  the  pool  side  edges. 
Further  study  is  needed  to  explain  this  discrepancy  between  theoretical  and 
experimental  results. 
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e.  Free  Surface  Deformation:  The  first  method  involves  solving  the  free 
surface  subroutine  in  conjunction  with  the  main  iteration  scheme.  By 
continuously  updating  the  free  surface  W  velocities  (and  consequently  the 
pressures)  deformation  of  the  surface  occurred.  A  higher  surface  tension 
value  (see  Table  3.3.1,  was  used  in  order  to  inhibit  the  free  surface 
deformation  from  jumping  over  the  phantom  nodes.  This  method  proved  to  be 
extremely  slow  due  to  the  convergence  requirements  of  the  surface  sub¬ 
routine.  A  much  quicker  method  used  to  visualize  the  surface  deformation 
resembles  a  perturbation  scheme.  The  converged  (flat  surface)  pressure  data 
was  used  as  initial  conditions  for  an  independent  surface  subroutine  calcu¬ 
lation.  The  new  surface  values  were  iterated  using  a  Poisson  pressure  formu¬ 
lation  until  conservation  of  mass  was  achieved  about  the  initial  surface 
position.  Instead  of  recalculating  the  free  surface  W  velocities  and 
continuing  the  main  iteration,  the  surface  values  were  simply  plotted,  as 
illustrated  in  Figure  3.3.22.  The  two  methods  yield  approximately  similar 
results  due  to  the  small  deformation  of  the  surface.  Method  2  is  more 
efficient  but  applies  to  small  surface  deformations  only.  It  should  be  noted 
that  the  "constant  volume"  convergence  criterion  was  used  at  the  melt  pool 
surface  because  exact  conditions  are  not  presently  known. 

Both  methods  show  the  following  trends.  Thermocapillarity  drives  the 
surface  fluid  radially  outward  at  extremely  high  velocities.  These  high 
velocities  displace  more  mass  from  the  central  (surface)  region  than  can  be 
replaced  by  the  recirculating  flow,  thus  causing  a  depression.  The  displaced 
mass  builds  up  at  the  solid/liquid  interface  causing  the  surface  to  bulge 
upward  where  it  is  then  forced  downward  into  the  molten  pool.  The  contact 


angle  of  the  liquid  metal  on  the  substrate  at  the  solid/liquid  interface  is 
unknown  and  currently  being  researched.  Dussan  [45}  indicates  that  the 
physics  of  contact  angle  are  very  poorly  understood  at  the  present  time. 

f.  Variable  Viscosity:  It  is  well  documented  that  some  material  proper¬ 
ties,  such  as  density,  remain  relatively  constant  while  others  such  as  surface 
tension  and  viscosity,  vary  considerably  with  temperature.  Variable  viscosity 
results  compared  with  constant  property  results  is  presented  in  Tables  3.34 
through  3.36  for  test  case  I,  Power  =  200  W,  =  0.010  m/s  which  is  illus¬ 
trated  in  Figs.  3.3.23  through  3.3.26.  (The  o  subscript  denotes  the  beam 
location  and  *  denotes  the  variable  viscosity  data.)  Temperature  and  velocity 
profiles  for  the  variable  viscosity  case  were  not  plotted  due  to  the  relative¬ 
ly  small  changes  in  the  velocities  and  pool  shapes  as  compared  to  the  constant 
property  case. 

An  exponential  curve  fit  was  used  to  model  the  viscosity  change  with 
temperature  in  the  melt  pool  region.  The  variation  of  viscosity  with  temper¬ 
ature  is  illustrated  in  Fig.  3.3.27.  Local  viscosity  varied  from  0.0055  kg/m- 
s  at  the  melt  temperature  to  approximately  0.0015  kg/m-s  near  the  concentrated 
heat  source.  Local  velocities  were  shown  to  increase  with  decreasing 
viscosity  (increasing  temperature),  the  most  dramatic  changes  taking  place 
near  the  free  surface  where  heat  is  being  convected  radially  outward  from  the 
concentrated  heat  source.  These  velocity  changes  were  minor  however, 
approximate  2-3  percent  on  the  average.  The  velocity  changes  were  smaller  in 
the  deeper  regions  of  the  melt  pool. 


g.  Electromagnetic  Forces:  When  a  current  is  passed  thru  the  scanning 
substrate  by  means  of  an  arc,  a  current  density  (J)  is  developed  in  the 
medium.  This  current  also  induces  a  magnetic  field  (B)  in  the  substrate  which 
interacts  with  the  J  field  to  form  electromagnetic  forces.  The  magnitude  of 
the  forces  is  determined  by  the  cross-product  of  the  J  and  B  terms  (JxB).  To 
illustrate  the  effects  of  the  JxB  forces  the  thermocapillary  influence  (for 
case  I)  was  removed  and  a  current  of  300  Amps  was  implemented  as  the  driving 
mechanism  of  the  flow.  As  seen  from  Figs.  3.3.28  through  3.3.31  the  JxB 
forces  drive  the  flow  radially  inward  and  downward  at  the  free  surface.  This 
downward  flow  impinges  upon  the  pool  bottom  and  recirculates  upward  along  the 
solid/liquid  interface,  this  recirculating  pattern  being  exactly  opposite  to 
that  of  the  thermocapillary  recirculation  pattern. 

The  JxB  forces  were  also  evaluated  along  with  the  surface  tension  driven 
forces  for  case  I.  A  current  of  100  Amps  was  used.  (See  Tables  3.37  through 
3.38).  (Once  again,  the  o  subscript  denotes  the  concentrated  heat  source 
center  and  the  *  denotes  the  variation  from  the  constant  property  case,  or  in 
other  words,  the  JxB  plus  surface  tension  data.)  Results  show  that  the 
radially  driven  surface  velocities  were  decreased  by  approximately  10  percent 
on  the  average  due  to  the  JxB  influence.  This  retardation  of  the  flow  tends 
to  shrink  the  pool  (radially)  and  deepen  the  melt  pool  bottom.  For  large 
enough  currents  this  inward  JxB  flow  will  predominate  over  the  outward  surface 
tension  flow  and  become  the  driving  mechanism  in  the  melt  pool. 

3.3.4  Conclusions 

From  the  numerical  solutions  the  following  trends  can  be  observed. 
Application  of  a  laser  onto  a  scanning  substrate  forms  a  melt  pool  which  has  a 


characteristic  length-to-width  aspect  ratio  of  approximately  1  and  a  charac¬ 
teristic  length-to-depth  aspect  ratio  of  approximately  7-12.  The  melt  pool 
isotherms  are  approximately  circular  in  the  (X-Y)  region  of  the  beam  center 
and  elongate  into  elliptic  shapes  away  from  the  beam.  The  elliptic  form  is 
due  to  the  influence  of  the  substrate  scanning  velocity  on  the  recirculating 
flow  field.  The  isotherms  tend  to  compress  near  the  solid/liquid  interface 
due  to  the  impinging  flow  while  the  spacing  of  the  isotherms  widens  near  the 
pool  bottom  due  to  the  reduced  velocities.  A  dimple  forms  at  the  melt  pool 
bottom  due  to  the  recirculating  nature  of  the  flow.  The  flow  is  driven  by 
surface  tension  effects  which  is  termed  thermocapillary  flow  due  to  the 
influence  of  temperature  on  surface  tension.  The  thermocapillary  velocities 
are  approximately  two  orders  of  magnitude  greater  the  substrate  scanning 
velocity  which  signifies  the  importance  of  surface  tension.  These  thermo¬ 
capillary  velocities  (on  the  order  of  1  m/s)  are  a  direct  consequence  of  the 
severe  temperature  gradients  imposed  on  the  substrate  by  the  laser.  These 
temperature  gradients  are  directly  proportional  to  the  cooling  rates 
surrounding  the  melt  pool  which  influences  the  microstructure  of  the 
resolidified  substrate.  Preliminary  free  surface  calculations  show  the 
development  of  a  depressed  region  in  the  melt  pool  center  causing  mass  to 
bulge  upward  at  the  solid/liquid  interface  before  being  recirculated  downward 
into  the  melt  pool.  Local  velocities  were  seen  to  increase  with  decreasing 
viscosity.  The  changes  were  small,  however,  approximately  2-3  percent  on  the 
average.  JxB  forces  induced  by  current  from  an  arc  tend  to  recirculate  the 
flow  opposite  to  that  of  the  surface  tension  flow,  and  may  predominate  if  the 
current  is  large  enough. 


These  results  give  important  insights  into  the  physical  phenomenon 
involved  in  laser,  arc  and  electron  beam  welding/alloying  of  metals.  Further 
investigation  of  this  subject  would  provide  information  on  controlling  the 
microstructure  of  the  weld  heat-affected-zone  (HAZ)  by  specifying  laser 
parameters  as  well  as  yielding  further  physical  insights  such  as  contact  angle 
and  oxidation  effects. 

Suggestions  for  program  improvement  include:  1)  the  incorporation  of  a 
variable  grid  scheme  for  finer  resolution  inside  the  melt  pool  boundary  layers 
and  a  coarser  resolution  in  the  moving  substrate  temperature  solution  domain, 
2)  the  incorporation  of  an  acceleration  scheme  for  quicker  residual  con¬ 
vergence,  3)  the  advancement  of  the  free  surface  routine  which  would  include 
steep  surface  curvature  and  multiple  node  jumping  capabilities  and,  4)  melt 
pool  coordinate  transformation  for  easier  computational  analysis  in  con¬ 
junction  with  a  cylindrical  coordinate  system. 
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3.3.6  Melt  Pool  Isotherms  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  240  W,  U.  =  0.012  m/s) 


Figure  3.3.7  Surface  Temperature  Contours 

(Power  =  240  W.  U  =  0.012  m/s) 
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Figure  3.3.10  Melt  Pool  Isotherms  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  400  W.  \1  «  0.02  m/s) 
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Figure  3.3 


.14  Melt  Pool  Velocity  Field  in  the  Y-Z  Plane 
with  X  Located  at  the  Beam  Center 
(Power  =  240  W,  U.  =  0.012  m/s) 
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Figure  3.3.15  Melt  Pool  Velocity  Field  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  240  W,  =  0.012  m/s) 
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Figure  3.3.16  Melt  Pool  Velocity  Field  in  the  X-Z  Plane 
with  Y  Located  at  the  Symmetry  Plane 
(Power  =  400  W,  Uj  =  0.02  m/s) 
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Figure  3.3.20  Surface  Temperature  versus  Dlstance(X)  for 

l.ocations(Y)  Away  from  the  Beam  Center 
(Power  »  400  W,  *  0.02  m/s) 
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Figure  3.3.21  Symmetry  Plane  Temperature  versus  Distance(X)  for 
Various  Locations(Z)  Beneath  the  Domain  Surface 
(Power  =  400  W,  U.  =  0.02  m/s) 
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Figure  3.3.22  Restricted  Free  Surface  Deformation  Based 
on  Local  Pressure  Only 
(Power  *  240  W,  U_  =  0.012  m/s) 


Figure  3.3.24  Melt  Pool  Isotherms  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  200  W.  Uj  =  0.010  m/s) 
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Figure  3.3.25  Melt  Pool  Velocity  Field  in  the  X-Z  Plane 
with  Y  Located  at  the  Symmetry  Plane 
(Power  =  200  W.  U  =  0.010  m/s) 


jc  Figure  3.3.26  Melt  Pool  Velocity  Field  in  the  X-Y  Plane 

with  Z  Located  at  the  Material  Surface 
(Power  =  200  W,  U5  =  0.010  m/s) 
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Figure  3.3.27  Variation  of  Viscosity  with  Temperature 


Figure  3.3.28  Melt  Pool  Isotherms  1n  the  X-Z  Plane 
with  Y  Located  at  the  Symmetry  Plane 
(Power  =  200  W,  =  0.010  m/s) 

(JxB  Forces  Only) 
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Figure  3.3.30  Melt  Pool  Velocity  Field  in  the  X-Z  Plane 
with  Y  Located  at  the  Symmetry  Plane 
(Power  =«  200  W,  U-  =  0.010  m/s) 

(JxB  Force  Only) 
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Figure  3.3.31  Melt  Pool  Velocity  Field  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  200  W,  =  0.010  m/s) 

(JxB  Forces  Only) 


Table  3.3.4  U-Velocity  Variations  Due  to  Variable  Viscosity 
Along  the  Surface  of  the  X-Z  Plane  with  Y  Located 
at  the  Symmetry  Plane  (Power  =  200  W,  =  0.010  m/s) 


(X-Xo) 

U 

U* 

aU 

% 

(mm) 

(m/s) 

(m/s) 

(m/s) 

-.320 

.000 

.000 

.000 

.000 

-.280 

-1.283 

-1.304 

.021 

1.624 

-.240 

-2.566 

-2.608 

.042 

1.625 

-.200 

-1.947 

-1.983 

.036 

1.867 

-.160 

-1.327 

-1.358 

.031 

2.335 

-.120 

-1.175 

-1.192 

.017 

1.467 

-.080 

-1.023 

-1.026 

.003 

.341 

-.040 

-.489 

-.497 

.008 

1.595 

.000 

.044 

.032 

-.012 

~n.S22 

.040 

.371 

.377 

.006 

1.691 

.080 

.697 

.722 

.025 

3.536 

.120 

1.048 

1.070 

.022 

2.110 

.160 

1.399 

1.418 

.020 

1.398 

.200 

1.486 

1.518 

.032 

2.159 

.240 

1.573 

1.618 

.045 

2.835 

.280 

2.255 

2.217 

-.038 

-1.679 

.320 

2.936 

2.816 

-.120 

-4.097 

.360 

1.468 

1.408 

-.060 

-4.098 

.400 

.000 

.000 

.000 

.000 

Table  3.3.5  V-Velocity  Variations  Due  to  Variable  Viscosity 
Along  the  Surface  of  th  Y-Z  Plane  with  X  Located 
at  the  Beam  Center  (Power  =  200  W,  =  0.010  m/s) 


(Y-Yo) 

V 

V* 

av 

% 

SI 

(mm) 

(m/s) 

(m/s) 

(m/s) 

-.320 

.000 

.000 

.000 

.000 

-.280 

-1.299 

-1.337 

.038 

2.916 

-.240 

-2.599 

-2.674 

.076 

2.916 

-.200 

-1.955 

-1.973 

.018 

.930 

-.160 

-1.311 

-1.271 

-.039 

-3.008 

-.120 

-1.102 

-1.072 

-.030 

-2.719 

«■ 

-.080 

-.893 

-.872 

-.020 

-2.296 

-.040 

-.447 

-.436 

-.011 

-2.487 

.000 

.000 

.000 

.000 

.000 

•  •. 

.040 

.447 

.436 

-.011 

-2.487 

.080 

.893 

.872 

-.020 

-2.296 

k 

.120 

1.102 

1.072 

-.030 

2.719 

.160 

1.311 

1.271 

-.039 

-3.008 

.200 

1.955 

1.973 

.018 

.930 

.240 

2.599 

2.674 

.076 

2.916 

.280 

1.299 

1.337 

.038 

2.916 

.320 

.000 

.000 

.000 

.000 

Table  3.3.8  Variation  of  V-Velocities  Due  to  Electromagnetic 
Forces  (JxD  with  Surface  Tension)  in  the  X-Y  Plane 
with  Z  Located  at  the  Material  Surface 
(Power  =  200  W,  =  0.010  m/s.  Current  =  100  Amps) 
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3.4  Experimental  Verification 
3.4.1  Introduction 

In  this  section,  two  experiments  are  described.  These  experiments 
are  performed  to  examine  the  validity  of  the  models.  It  should  be  kept  in 
mind  that  there  is  no  thermophysical  properties  (not  to  mention  the  temper¬ 
ature  dependence)  data  available  for  high  or  moderately  high  melting  point 
materials  such  as  steel  and  aluminum.  In  view  of  this,  exact  match  of  experi¬ 
mental  and  theoretical  predictions  are  unlikely.  However,  trends  of  experi¬ 
mental  and  theoretical  predictions  should  agree  to  each  other.  The  first 
experiment  is  to  compare  the  aspect  ratio  of  the  pool  shape  (width/depth)  of 
the  experimental  and  theoretical  predictions.  It  can  be  seen  from 
Eq.  (3.1.19)  that  both  Reynolds  number  and  Marangoni  number  are  directly 
proportional  to  the  power  absorbed.  These  parameters  therefore  serve  the 
purpose  of  the  comparison.  The  second  experiment  is  to  obtain  the  magnitude 
of  the  recirculating  flow  and  the  cooling  rate.  This  is  down  by  x-ray  shadow 
graph  using  a  real-time  radiography  system. 

3.4.2  Experimental  Procedure 

a.  Micrograph  of  Resolidified  Region;  A  schematic  diagram  of  the 
experimental  set  up  is  shown  in  Fig.  3.4.1.  A  10  kW  CO2  continuous  laser  is 
used.  The  materials  used  are  AISI  1016  mild  steel.  The  steel  bar  is  of 
dimension  150  mm  long,  50  mm  wide,  and  6  mm  thick.  The  surface  being 
irradiated  by  the  laser  is  ground.  The  surface  is  cleaned  by  using  alcohol 
right  before  the  run  to  avoid  any  unnecessary  containment. 

Samples  for  the  cross-section  of  laser  melt  were  cut,  hand  polished,  and 
etched  for  optical  microscopy.  The  melted  region  is  identified  and  the  aspect 
ratio  of  the  melted  region  is  estimated. 


b.  X-ray  Shadow  Graph:  A  schematic  diagram  of  the  experimental  set  up  is 
shown  in  Fig.  3.4.2.  An  arc  is  used  as  the  heat  source.  The  materials  used 
are  aluminum.  Tungsten  particles  are  used  as  tracers  to  provide  contrast. 
Microfocus  x-ray  is  used  to  produce  a  focused  x-ray  beam.  This  x-ray  beam  is 
directed  through  the  test  section.  An  image  intensifier  is  located  on  the 
back  to  receive  the  image.  This  image  is  processed  and  stored  on  video 
tape.  Because  of  density  difference  between  aluminum  and  tungsten,  the 
tungsten  particles  can  be  traced  as  they  travel  with  the  motion  of  the  molten 
aluminum.  Density  difference  due  to  temperature  variation  can  also  be 
recorded. 

3.4.3  Results  and  Discussion 

a.  Micrograph  of  Resolidified  Region:  The  experimental  runs  for 
various  process  parameters  -  laser  power,  beam  diameter,  and  scanning  speed  - 
are  tabulated  in  Table  3.4.1.  A  picture  of  the  microstructure  of  case  #6  is 
shown  in  Fig.  3.4.3.  The  melted  region  can  be  distinguished  from  the  unmelted 
region.  The  aspect  ratios  based  on  these  micrographs  are  calculated  and 
tabulated  in  Table  3.4.1.  It  can  be  seen  that  the  aspect  ratio  increases  with 
respect  to  laser  power  which  is  consistent  with  the  theoretical  prediction. 
The  experimental  aspect  ratio  varies  from  3.375  to  4.0.  The  theoretical 
models  predicts  4.0  to  5.0  for  the  aspect  ratio. 

b.  X-ray  Shadow  Graph:  In  Fig.  3.4.4,  two  successive  frames  of  the  x-ray 
shadow  graph  are  shown.  The  locations  of  the  center  line  and  tungsten 
particle  are  shown  by  the  vertical  lines.  The  left  line  directly  underneath 


the  arc  is  the  center  line.  The  right  line  is  the  location  of  the  tungsten 
particle.  Assuming  that  the  particle  travel  with  the  motion  of  the  molten 
materials,  the  speed  of  the  recirculating  flow  can  be  estimated.  It  is  found 
to  be  of  the  order  1  m/sec  which  is  of  the  same  order  of  magnitude  o'"  the 
theoretical  prediction. 

A  sequence  of  four  frames  at  0,  1,  5,  and  9  seconds  after  the  arc  is 
turned  off  are  restored  by  using  a  frame  grabber.  The  first  frame  is  used  as 
a  reference.  Three  consecutive  frames  each  four  seconds  apart  are  produce 
from  the  next  three  frames  and  plotted  in  Fig.  3.4.5  through  3.4.7.  Each 
frame  is  produced  by  subtracting  the  intensity  from  the  reference  frame.  The 
difference  in  intensity,  which  is  caused  by  different  density  and  hence 
different  temperature,  is  color  enhanced.  The  change  in  the  color  is  there¬ 
fore  a  measure  of  the  cooling  rate.  It  can  be  observed  that  the  cooling  rate 
is  non-uniform  within  the  molten  region.  It  has  a  higher  cooling  rate  on  the 
surface  than  near  the  bottom.  This  is  consistent  with  the  theoretical 
prediction. 

3.4.4  Conclusion 

Two  experiments  are  performed  to  examine  the  validity  of  the 
theoretical  models.  The  first  experiment  compares  the  aspect  ratio  of  the 
molten  pool.  It  is  found  that  the  theoretical  and  experimental  results  have 
the  same  trend.  An  increase  in  laser  power  increases  the  aspect  ratio.  The 
aspect  ratio  of  the  two  predictions  are  within  20  percent.  The  second 
experiment  is  to  obtain  a  measure  of  the  velocity  field  and  cooling  rate.  It 
is  found  that  the  recirculating  flow  is  of  the  order  1  m/sec.  The  cooling 
rate  within  the  molten  region  is  found  to  be  non-uniform.  It  is  higher  on  the 
surface  and  lower  near  the  bottom. 
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Figure  3.4.1  Schematic  Diagram  of  the  Experimental 
Set  Up  -  Resolidified  Region 
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Table  3.4.1  Numerical  Values  of  Process  Parameters  of  the  Experimental  Runs 


Beam  Radius  0.5  mm 


Case  # 

Scanning  Speed 

Laser  Power 

Aspect  Ratio 

1 

25.0 

mm/sec 

2  kW 

3.375 

2 

25.0 

mm/sec 

3  kW 

3.51 

3 

25.0 

mm/sec 

4  kW 

3.62 

4 

50.0 

mm/sec 

3  kW 

3.66 

5 

50.0 

mm/sec 

4  kW 

3.79 

6 

50.0 

mm/sec 

8  kW 

4.0 

127 


4.  SUMMARY  AND  CONCLUSIONS 

Modeling  efforts  carried  out  in  the  first  two  years  of  this  ONR  program 
improved  the  general  understanding  of  the  process  and  revealed  the  different 
regimes  of  the  flow  as  shown  in  Fig.  1.1.  It  also  drew  our  attention  to  the 
singularities  and  limiting  behavior  of  the  flow  which  will  be  addressed  in  our 
future  work.  The  important  observation  and  correlations  focus  the  present 
study  can  be  summarized  as  follows. 

1.  Perturbation  Model 

The  physics  of  the  process  changes  in  the  presence  of  thermocapillary 
convection.  It  becomes  convection  dominated.  This  change  the  pool  geometry 
dramatically  (150  percent  increase  in  the  aspect  ratio  (width/depth)  as 
compared  to  the  conduction  case).  The  cooling  rate  along  the  solidifying 
surface  is  found  to  be  non-uniform.  This  in  turn  give  rise  to  the  non-uniform 
resolidified  microstructure.  The  recirculating  flow  in  the  molten  region 
gives  a  possible  explanation  of  solute  redistribution. 

2.  Stagnation  Flow  Analysis 

This  model  represents  the  dominant  feature  in  the  central  region  of  the 
melt  pools  in  various  materials  processing  processes.  This  model  clarifies 
the  scaling  laws  for  thermocapillary  convection  in  a  deep  melt  pool.  Single 
and  concise  expressions  are  obtained  for  the  physical  quantities  such  as 
velocity,  temperature,  and  viscous  and  thermal  boundary  layers.  Prandtl 
number  dependence  is  also  characterized.  Good  agreements  are  obtained  in 
comparison  with  the  basic  flow  of  the  perturbation  model.  The  simple  and 
concise  expressions  can  therefore  be  used  to  estimate  velocity  and 
temperature. 


This  model  provides  the  full  three-dimensional  realistic  calculations  of 
the  process.  Results  of  both  thermocapillary  and  Lorentz  flow  are  presented 
and  discussed.  It  is  found  that  the  thermocapillary  flow  tends  to  shear  the 
liquid  metal  to  the  pool  edge  while  Lorentz  force  causes  the  flow  in  the 
opposite  direction.  Depending  on  their  magnitudes  one  might  dominate  over  the 
other,  preliminary  results  of  the  free  surface  shape  due  to  fluid  motion  is 
also  presented  and  discussed. 

4.  Experimental  Verification 

Two  experiments  were  performed  to  examine  the  validity  of  the  models.  It 
is  found  that  both  theoretical  models  and  experiments  predict  the  same  trend 
in  aspect  ratio  (width/depth).  The  aspect  ratio  increases  with  laser  power. 
It  is  also  found  that  the  recirculating  flow  within  the  molten  pool  is  of  the 
order  1  m/sec  which  coincides  with  theoretical  prediction.  The  cooling  rate 
without  the  molten  region  is  found  to  be  non-uniform.  It  is  higher  on  the 
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